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ABSTRACT 



Analysis of INTEGRAL Core Program and public Open 
Time observations has recently provided a sample of 60 
extragalactic sources selected in the 20-100 keV band 
above a flux of 1.5 10~ n erg cm~ 2 s _1 . As this band 
probes heavily obscured regions/objects, i.e. those that 
could be missed in optical, UV, and even X-ray surveys, 
our sample offers the opportunity to study the extragalac- 
tic sky from a different point of view with respect to sur- 
veys at lower energies. We present an update of our anal- 
ysis, including the first sample of AGNs detected above 
100 keV. We also discuss the results of follow up observa- 
tions performed at optical and X-ray frequencies with the 
aim of classifying our objects and studying the effects of 
intrinsic absorption in gamma-ray selected AGNs. The 
average redshift of our sample is 0.134 while the mean 
20-100 keV luminosity in Log is 43.84; if blazars are 
excluded these numbers become 0.022 and 43.48 respec- 
tively. Defining an absorbed object as one with Nu above 
10 22 atoms cm~ 2 , we find that absorption is present in 
60% of the objects with at most 14% of the total sample 
due to Compton thick active galaxies. Almost all Seyfert 
2s in our sample are absorbed as are 24% of Seyfert 
Is. We also present broad-band spectral information on 
a sub-sample of the brightest objects: our observations 
indicate a mean photon index of T=1.8 spanning from 
30-50 keV to greater than 200 keV. Finally, we discuss 
the LogN/LogS distribution in the 20-100 and 100-150 
keV bands derived from our sample. The present data 
highlight the capability of INTEGRAL to probe the ex- 
tragalactic gamma-ray sky, to discover new AGNs and to 
find absorbed objects. 

Key words: galaxies: active - gamma rays: observations: 
surveys. 



1. INTRODUCTION 



Quantifying the fraction of AGNs missed by surveys 
which are affected by selection due to absorption is nec- 
essary if we want to fully understand the accretion his- 
tory of the Universe and study a population of objects 
(absorbed AGNs) so far poorly explored. Furthermore, 



the measurement of the primary AGN continuum and 
its cut-off energy is crucial for understanding emission 
models and discriminating between them. Both informa- 
tion (spectral shape of AGN and column density distribu- 
tion) are key parameters for estimating the contribution 
of AGNs to the X-ray cosmic diffuse background and for 
testing current unified theories. Only high energy obser- 
vations (above 10 keV) can provide an unbiased knowl- 
edge of the column densities and high energy cut-off dis- 
tributions of AGNs in the local Universe. 
A step forward in both these issues is now provided by 
INTEGRAL/IBIS and SWIFT/BAT which are surveying 
a great fraction of the sky above 20 keV with a sensitiv- 
ity better than a few mCrab and a point source location 
accuracy of the order of 1-3 arcminutes depending on the 
source strength. These two surveys are complementary to 
each other, as the first one concentrates mostly on map- 
ping the galactic plane, while the second covers the high 
galactic latitude sky so that together they will provide the 
best yet sample of AGNs selected in the gamma-ray band. 
We have recently presented a catalogue of ~ 60 IBIS 
AGNs selected in the 20-100 keV band above a flux 
limit of 1.5xl0 -11 erg cm~ 2 s _1 (|3] or paper I). The 
analysis was performed using INTEGRAL Core Program 
and public Open Time observation performed up to April 
2005. Most sources in the sample were Seyfert galaxies 
almost equally divided between type 1 and 2 objects, 6 
were blazars but 14 were still unclassified at the time of 
this first publication. Despite an initial attempt to study 
the role of absorption in the sample, a significant fraction 
(about 30%) of the objects did not have archival X-ray 
spectra to allow a full assessment of the column density 
distribution. Furthermore due to the highly inhomoge- 
neous coverage of the INTEGRAL/IBIS survey a deter- 
mination of the LogN-LogS relation had to await a proper 
study of the various corrections to be applied. 
Here, we present an update on the first 20-100 keV 
IBIS catalogue of AGNs, with a number of identi- 
fication/classification being provided in the meantime 
through optical spectroscopy; on the basis of further work 
we have eliminated a few sources from the initial list and 
added a couple of new objects. We also use available X- 
ray data to obtain information on the column density of a 
number of AGNs in the catalogue. Given the higher num- 
ber of new information available, we have re-estimated 
the main conclusions of 01 ■ Furthermore on the basis of 
this update, we have evaluated the LogN-LogS relation 



in the 20-100 keV band. In the meantime, we have also 
performed broad band spectral analysis of a subsample 
of objects observed in the first INTEGRAL extragalactic 
survey |2] using IBIS data in combination with archival 
X-ray data from other missions and have provided infor- 
mation on their primary continuum shape including the 
cut-off energy. Finally, we report the first catalogue of 
AGNs detected above 100 keV by INTEGRAL and pro- 
vide the first extragalactic number counts in the 100-150 
keV band. 



2. CATALOGUE UPDATE AND OPTICAL CLAS- 
SIFICATION 



Of all the excesses reported in paper I, 3 objects re- 
sulted a posteriori to be artifacts of the imaging pro- 
cesses mainly due to projection effects; these 3 sources 
(IGRJ 13000+2529, IGR J 13057+2036 and UGC3142) 
have been removed from the initial sample. Two ob- 
jects which were not listed but were detected in the sur- 
vey, were found to be AGNs a posteriori through opti- 
cal follow up observations ( IVm and references therein); 
due to the limited information available on these two ex- 
cesses (IGR J17488-3253 and IGR J17513-2011) their 
AGNs nature was difficult to assess prior to optical spec- 
troscopic work. It is important to stress that some other 
unidentified sources may turn out to be active galaxies 
but due to their location near the galactic plane it is dif- 
ficult to assess their nature without a proper optical fol- 
low up work. We have estimated for these two objects 
the total exposure and 20-100 kev flux and luminosity as 
done in paper I for the other sources; ISGRI coordinates 
can instead be found in |7]. We obtained an exposure of 
737.0 and 571.8 ks, a flux of 2.18±0.17 and 2.04±0.22 
mCrab and a luminosity of 3 and 15.5 x 10 43 erg cm -2 
s- 1 for IGR J17488-3253 and IGR J17513-201 1 respec- 
tively. We remind that between 20-100 keV, 1 mCrab 
corresponds, for a Crab like spectrum, to 1.6 x 10~ n erg 
cm _2 s _1 which is a value close to our detection limit; the 
fluxes have been converted to gamma-ray luminosities as- 
suming Ho=71 Km sec -1 Mpc -1 and qo=0 |21]. Finally 
one source (IGR J 16 194 -2810) has been removed as it 
is not clear at the present stage from our own follow up 
measurements (Masetti, private communication) if it is a 
galactic or an extragalactic source. A number of objects 
in the initial catalogue have been optically classified and 
confirmed to be active galaxies; at the moment only 5 
objects have no classification nor redshift available. For 
all the objects in this initial survey, we report in Table 1 
an update of the source classification and redshift either 
available from the literature or from our own work (" 11711 
at this conference and references therein); in some cases 
the optical classification has been revised allow interme- 
diate Seyfert types to be also sampled. We have also up- 
dated information on the column density by using avail- 
able X-ray data. In particular we have made use of public 
SWIFT/XRT observations and performed our own spec- 
tral analysis in order to estimate photon index and column 
density for a set of 8 sources. More detailed results of this 
analysis are presented in the contribution to this confer- 



10 ; 




0.001 0.01 0.1 1.0 10. 

Redshift 



Figure 1. 20-100 keV Luminosity versus redshift for opti- 
cally classified AGNs in the present survey; straight line 
corresponds to the IBIS survey limit 

ence by Landi et al. [ 12]. We have also used XMM data 
obtained very recently in conjunction with INTEGRAL 
data to quickly estimate the source column density. A 
more in depth analysis of these objects is on going but 
since we are only interested in estimating if the column 
density is above or below 10 22 at cm~ 2 a simple mod- 
eling of the data (power law absorbed both by galactic 
and where required by intrinsic absorption) is sufficient 
to provide this information. Only in the case of IC4518A 
we find that an intrinsic column density above 10 22 at 
cm~ 2 is required by the data. Results becoming avail- 
able in the literature have also been used to update our 
database 12311 . Only 9 objects have no column density 
estimate at this stage. 



3. RESULTS UPDATE 



Our extragalactic survey comprises after this revision 60 
objects securely identified with AGNs. For those objects 
with known distance, we plot in figure 1 the gamma-ray 
luminosity against redshift, to show the large range in 
these parameters sampled by the present survey. ^From 
this figure it is also evident that our sensitivity limit is 
around 1.5 x 10~ n erg cm _2 s _1 (straight line in the 
figure). Within the overall sample, 49 objects are clas- 
sified as Seyfert galaxies, 6 are blazars and only 5 (8% 
of the sample) are still unclassified. Within the sample of 
Seyfert galaxies, 26 objects are of type 1-1.5 while 23 are 
of type 2, i.e. a ratio 1:1, which illustrates the power of 
gamma-ray surveys to find narrow line AGN. About 10% 
of our sample is made of radio loud objects. The average 
redshift of our sample is 0.134 while the mean 20-100 
keV luminosity in Log is 43.84; if blazars are excluded, 



Table 1. AGNs sample 
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Tab 1 - Continued 



Source 


Type 


z 


Log Nh 


IGR J 16 119-6036 


SI 


0.016 


- 


IC4518-A 


S2 


0.016 


23.3$ 


NGC6221 6 


S2 


0.005 


22.0 


IGR J16558-5203 


S1.2 


0.054 


<20.0 


IGR J 17204-3554 


- 




23.1 


IGR J 18249-3243 


- 




- 


IGR J20 187+4041 


- 




23.3 


IGR J20286+2544 


S2 


0.014 


23.6-24.0 


IGR J2 1178+5139 









f Landi et al. these proceedings 

X Panessa et al. in preparation and De Rosa et al. in 
preparation 

★ Piconcelli et al. (2006), also absorption partially 
covering the source 



these numbers become 0.022 and 43.48 respectively. 
The column density distribution for our sources is pre- 
sented in figure 2. Assuming 10 22 atoms cm~ 2 to be the 
dividing line between absorbed and unabsorbed sources, 
we find that absorption is present in nearly 60% of the 
sample of objects having a column density estimate. This 
result is in line with the Swift findings 1 16] and also with 
other INTEGRAL AGNs surveys jj, j2j! 
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Figure 2. Column density versus 20-100 keV Luminosity 
for AGNs in the survey with known intrinsic absorption 

Within the present sample we find that 6 possibly 7 ob- 
jects are Compton thick. NGC4945, Circinus galaxy, 
NGC1068 and MKN3 are well known examples of this 
class |H. To this set we add ESO 138-G01 1 d, 

1 analysis of various energy band mosaics indicates detection of both 
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Figure 3. Contours of column densities in 4U 1344-60 (one fully versus one partially covering the source) and in MKN 6 
( in this case both partially covering the source ) ASCA and BeppoSAX data clearly indicating Nh variability. 



possibly NGC6300 since this is an object which moves 
from thick to thin in time IU9I1 and very likely IGR 
J20286+2544 (H, O- Therefore we estimate that at 
most 14% of the sample is Compton thick. Within the 
sub-sample of 22 Seyfert 2s with known Nh, we find that 
95 % are absorbed and 30% are Compton thick; this is in 
line with previous estimates of the column density dis- 
tribution of type 2 objects based on X-ray data l24ll . (jj- 
Interestingly we also find that 24% of type 1-1.5 objects 
with Nu values are absorbed. In some of these absorbed 
type 1 objects, the absorption is complex with multiple 
column density obscuring partially and/or totally the cen- 
tral nucleus; in same cases these absorptions could also 
vary in time. Typical examples of this complexity are 
MKN 6, NGC4151 and 4U 1344-60 (see |H and refer- 
ences therein). In figure 3 we show the contours of the 
column densities for 4U 1344-60 and MKN 6 P . 
Similar sources are also found in the SWIFT/BAT survey 

In figure 4 we show the column density of our sources as 
a function of IBIS luminosity. Although from this figure 
there is little evidence of any strong correlation, the few 
high luminosity objects in the survey all have low absorp- 
tion. However, since high luminosity objects are mostly 
blazars, this evidence is not conclusive and more data are 
required to firmly establish this observational evidence. 
More in depth studies of the present sample require op- 
tical classification of all objects and a detailed analysis 
of their broad band behavior (particularly in the X-ray to 
gamma-ray band) in order to fully understand the role of 
absorption : in particular we want to establish if there is 
any other Compton thick object in the sample and how is 
the absorption (complex or not) in Seyfert 1-1.5 galaxies. 
Furthermore when all sources are properly characterized 
it will be possible to estimate the luminosity function and 
the role of gamma-ray selected AGNs with respect to the 
cosmic diffuse background (see also l25ln . 
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Figure 4. Column density versus 20-100 keV luminosity 
for AGNs in the present sample 



4. STUDY OF THE PRIMARY CONTINUUM 



The high energy emission of AGNs is often to first 
order well described by a power law of photon index 
1.8-2.0, extending from a few keV to over 100 keV; 
at higher energies there is evidence of an exponential 
cut-off. Secondary features, which are also commonly 
present, are considered to be the effects of reprocessing 
of this primary continuum and are relatively well un- 
derstood. Modeling of high energy AGN spectra has so 
far generally focused on how to reproduce and explain 
the observed primary continuum shape: however, while 
the photon index distribution has been well investigated 
fT8Tl . observational results on the cut-off energy have 



Figure 5. High Energy cut-off vs photon index for MCG-5-23-16 (right) and GRS 1734-292 (left). 



so far been limited by the scarcity of measurements 
above 10-20 keV, with most information coming from 
BeppoSAX broadband spectra. Analysis of type 1 and 2 
AGN 1 22] provides evidence for a wide range of values 
in the cut-off energy, spanning from 30 to 300 keV and 
further suggests a possible trend of increasing cut-off 
energy as the power law index increases; it is not clear 
however if this effect is due to limitations in the spectral 
analysis or if it is intrinsic to the source population 
sampled. 
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Figure 6. Power law photon index versus cut-off energy 
for the sample of AGNs analyzed by Soldi et al. 2005 and 
Molina et al 2006. 

Direct INTEGRAL measurements enable us to obtain 
further observational results on the primary continuum 
and high energy cut-off thus providing more refined pa- 
rameters for AGN modeling. The sample of 12 AGNs 
first detected by INTEGRAL H can be taken as a case 



study: they are representative of the larger sample of 
AGNs presented here and are detected with adequate 
statistics to allow more in depth analysis. In particular, 
they are ideal objects for the study of the primary power 
law component, the presence of any high energy cut-off 
and the existence of a relation (if any) between these two 
parameters. Spectral analysis has recently been reported 
for this entire sub-sample using INTEGRAL data alone 
or in conjunction with data from other satellites 12611 . 
JToh . |H. Excluding PKS1830-211 which is a blazar, 
this sub-sample is primarily made of Seyfert galaxies. 
Indeed PKS 1830-211 shows a flat spectrum and a cut- 
off likely located at MeV energies as typical of a low 
frequency peaked or red blazar llOfl . For the remaining 
sources, the mean photon index is 1.83±0.07, while the 
cut-off energy ranges from 30-50 keV to greater than 200 
keV. In most cases the cut-off energy is well constrained 
at or below 100 keV: the contour plot of photon index 
versus cut-off energy for MCG-05-23-16 and GRS 1734- 
292 are shown in figure 5 and indicates the ability by high 
energy observations of putting tight constraints on both 
parameters . Combining all our data, we can also test 
the correlation found by a number of authors between the 
photon index and the cut-off energy. In figure 6, we plot 
the photon index versus the high energy cut-off for those 
Seyfert galaxies in the sub-sample for which we have in- 
formation on these two parameters. It is clear from this 
figure that our data do not show the correlation found 
by 1 22]: the values for the photon index, in fact, cluster 
around 1.8, while the values of the cut-off energy range 
over a wide interval. We can conclude that the correlation 
between the high energy cut-off and the photon index re- 
mains to be proved and that a diversity in cut-off energy 
is most likely a property of Seyfert galaxies. INTEGRAL 
will keep observing these AGNs and new ones for the du- 
ration of the mission and so we expect that more spectral 
data on this class of objects will become available; this 
will provide further information on the high energy cut- 
off in AGNs and a deeper insight into the primary contin- 
uum emission. 
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Figure 7. HR1 (40-100 keV/20-40 keV) versus HR2 (100- 
150 keV/40-100 keV)for the lOAGNs detected above 100 
keV 



5. AGNS DETECTED ABOVE 100 KEV 



Recently Bazzano et al. |4j] have presented the first cen- 
sus of IBIS objects with detection above 100 keV. This 
high energy survey includes 10 AGNs with detection 
above 4 sigma level in the 100-150 keV band of which 
only one (Cen A) is also visible in the 150-300 keV band. 
All objects are also detected at energies below 100 keV. 
Of these 10 AGNs, 2 are blazars and 8 are Seyfert galax- 
ies; however this time the ratio of type 1 versus type 2 
objects is 2:6, with Seyfert 2 outnumbering Seyfert 1 by 
a factor of 3, which is interestingly close to values found 
in optical spectroscopic surveys (3, full . However, 
this could be simply due to the limited number of AGNs 
present in this first extragalactic survey above 100 keV. In 
Figure 7, we plot for these 10 AGNs the hardness ratios 
HI versus H2, where HI is defined as the 40-100 keV to 
20-40 keV flux and HR2 as the ratio of the 100-150 keV 
to the 40-100 keV flux. It is clear from this figure that 
the 10 AGNs detected at high energies have photon in- 
dex > 1.5 with some clearly showing a softening of the 
spectrum probably due to the presence of a cut-off at or 
around 100 keV as discussed in the previous section. 



6. AGN LOGN-LOGS IN THE 20-100 KEV AND 
100-150 KEV BANDS 



In this section, we estimate the number counts relation for 
AGNs in two bands : 20-100 keV and 100-150 keV. Even 
though the number of sources in the second band is low, a 
LogN-LogS estimate is nevertheless a useful information 
to have if one wants to estimate the total number of ob- 
jects which may become visible as the survey progresses. 
However before doing any evaluations, there are a num- 
ber of complications involved in forming this relationship 
that must be taken into consideration. First, the sensitiv- 
ity limit in any direction depends not only on the expo- 
sure at that point, which is extremely non-uniform due to 
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Figure 8. The la background variation as a function of 
exposure in the 20-100 keV band. 



the pointing strategy of INTEGRAL 01, but also on sys- 
tematic variations/structures in the background. Another 
problem is the fact that our survey is not truly serendipi- 
tous as some of our objects have been the targets of IN- 
TEGRAL dedicated observations. Finally, our sample is 
probably incomplete again due to the peculiar nature of 
INTEGRAL observing strategy. Given these caveats, we 
first convert exposure to 1 a limiting sensitivity the flux 
errors and effective exposures for the sources found in 
each band, (see fig 8 for the case where the errors refer to 
20-100 keV fluxes). From this we obtain the sky cover- 
age as a function of 5a (4a) 20-100 keV (100-150 keV) 
limiting flux as shown in figure 9, again for the case of 
20-100keVband. 

We then create the 'raw' logN-logS for the 60 AGNs 
detected in the 20-100 keV band (figure 10) and the 10 
found in the 100-150 keV band (figure 11). This distri- 
bution is then corrected for the exposure limit and sky 
coverage assuming a-priori an isotropic underlying dis- 
tribution. Both the 'raw' data and the corrected values 
normalized to full sky coverage are also shown in figures 
10 and 1 1 where the horizontal bars indicate the errors on 
the flux of the N th source. The best-fit power law rela- 
tionship and the (la) statistical limits are indicated by the 
solid lines. Unsurprisingly, we find that the best-fit power 
law is similar to, though rather flatter than an isotropic 
distribution. A deeper analysis will be performed as a 
greater number of sources and, most importantly, more 
uniform sky coverage become available. 

Finally these number counts can be used to estimate the 
contribution of local AGNs to the cosmic diffuse back- 
ground. The total AGNs emissivity measured in both 
bands above 1 mCrab accounts for approximately 1-3% 
of the total background intensity; extrapolation of our 



100 




5(7 Minimum Flux (mCrab) 



Figure 9. The percentage of sky coverage as a function of 
5 a limiting flux. 



LogN-LogS by a factor of 100 towards lower flux will 
account for 15-20% of the extragalactic background. In 
short we are just start to unveil the source population re- 
sponsible for the hardX/gamma-ray background. 



7. FUTURE EXPECTATIONS 



Future work will concentrate on the optical classifica- 
tion of objects in the present sample and of those new 
AGNs that will be detected in the future. We also want 
to progress in estimating the column densities for the 8 
remaining sources in the 20-100 keV catalogue and pre- 
pare to obtain X-ray data on the new sources. In parallel, 
we expect to concentrate on detailed analysis of the spec- 
tral properties of all AGNs in the present sample in order 
to define the primary continuum (photon index and cut- 
off energy) using INTEGRAL data in combination with 
archival X-ray information and data obtained though pro- 
posal requests to the Chandra, XMM and Swift observa- 
tories. In figure 12 we plot the number of AGNSs de- 
tected in our own previous INTEGRAL surveys 01, 0|, 
J3| as a function of science windows analyzed. We are 
currently working at the third INTEGRAL survey which 
involves the analysis of around 24000 science window. 
As evident in the figure we expect around 100 AGNs 
(even more as new part of the extragalactic sky will be 
observed for the first time) to be detected in this coming 
survey thus doubling the number of sources for which 
follow up optical and X-ray observations will be needed, 
investigation of their spectral properties required and sta- 
tistical studies performed. The future looks good for ex- 
tragalactic studies with INTEGRAL and busy for the ob- 
servers. 
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Figure 10. 20-100 keV full sky number-flux relationship 
for our sample 
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Figure 11. 100-150 keV full sky number-flux relationship 
for the sample of 10 AGNs reported in [ffl. 
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Figure 12. The increase in AGN numbers as a function of 
accumulated science windows. When the 3 rd IBIS/ISGRI 
survey becomes available around 1 10 AGN should be vis- 
ible. 
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